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I. INTRODUCTION 



Despite several decades of scrutiny, the intrinsic shape of the nucleon, a fundamental building block of our world, 
is still not fully resolved [1, 2]. Although the nucleon is experimentally easily accessible in exclusive and inclusive 
scattering its spectroscopic quadrupole moment is zero since it has spin J = 1/2. However, this does not mean 
that the nucleon is spherically symmetric, since it can have an intrinsic deformation. Deformation is a common 
phenomenon in nuclear and atomic physics. Quantum mechanically, a multiphoton coincidence experiment could 
determine that a J — ground state of a diatomic molecule or nucleus has a deformed shape. However, usually in 
electromagnetic probes of microscopic systems, we are constrained to make measurements associated with one-photon 
exchange, corresponding to a matrix element of a one-body operator. In the case of a diatomic molecule, the one-body 
charge density of the J = state is spherically symmetric, and cannot reveal the deformation that is present in the 
system. For J > 1/2, however, when a nuclear or atomic system is well approximated by a deformed intrinsic state, 
it is still possible to observe its deformation using a one-body electromagnetic operator. In this case, the quadrupole 
moment is non-zero in the laboratory frame if the state is deformed. For collective rotation of the deformed intrinsic 
state [3], the relation between the spectroscopic quadrupole moment, Q, measured in the laboratory frame and the 
intrinsic quadrupole moment, Qo, in the body- fixed intrinsic frame is given by 

_ 3ir 2 -j(j + i) 

Q - (J+l)(2J + 3) g ° ' (1) 

where J is the total angular momentum of the system in the laboratory frame, K is the projection of J onto the z-axis 
of the body-fixed intrinsic frame, and we have considered the sub-state with azimuthal quantum number M = J. In 
the previous example of the J = diatomic molecule, although Qo 7^ 0, Eq. (1) yields Q — so that the deformation 
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is unobservable. Similarly, in the case of a nucleoli with J = 1/2, Q is zero although Q may not be [4, 5]. Since the A 
has J = 3/2 a non-zero quadrupolc moment Q can be measured [6]. The electric and Coulomb transition amplitudes 
E2 and C2 between the J = 1/2 nucleon and its J = 3/2 resonance have the same property of revealing the presence 
of deformation in the nucleon, the A, or both. Therefore, as in experiment, in this work we detect deformation by 
measuring E2 and C2. 

In recent years, we have seen tremendous progress in experimental measurements of the subdominant quadrupolc 
amplitudes [7, 8, 9, 10, 11, 12, 13, 14], yielding very accurate results particularly at low q 2 . For a recent review 
of the experimental situation see Ref. [15]. These accurate measurements have motivated several recent theoretical 
studies both in lattice QCD [16, 17, 18] and in chiral effective theories [19, 20]. On the lattice, hadron deformation 
can also be studied by investigating directly the charge distribution using density-density correlators [4, 21, 22], and 
it was shown that the rho meson is deformed [23]. The A shows deviations from a spherical shape albeit with large 
statistical errors. The issue of deformation of the A using density-density correlators is under study with improved 
lattice techniques [24, 25]. For concise reviews see [5, 26]. 

The focus of this work is a calculation of the nucleon - A electromagnetic transition form factors within QCD. We 
use two different types of simulations: The first uses two degenerate flavors of dynamical Wilson fermions and the 
second a hybrid action. The hybrid action uses dynamical staggered sea quarks, with two degenerate light quarks 
and one fixed to the mass of the strange quark. These dynamical quark configurations are produced by the MILC 
collaboration [27] and represent a state-of-the-art simulation of the QCD vacuum with the three lightest flavors of 
quarks taken into account. For the light valence quarks we use two degenerate domain wall fermions. This approach 
has been used successfully to evaluate nucleon structure functions [28, 29, 30, 31, 32, 33, 34] and the N to A axial- 
vector form factors [35, 36, 37, 38]. In this work, we compare results calculated using the hybrid action to the results 
obtained using two degenerate flavors of dynamical Wilson fermions. Given the different lattice artifacts involved in 
the two approaches, agreement between them provides a consistency check of our lattice methodology. Comparison 
between the dynamical results and our quenched results [16, 39] probes pion cloud contributions. 

The first lattice study of the N to A electromagnetic transition form factors [40] was carried out in the quenched 
approximation with limited statistics. Although the mean value was negative, the statistical errors on the suppressed 
quadrupole amplitudes were large and a zero value could not be excluded. Nonetheless, this pioneering work set up 
the methodology for a more elaborate study that would soon become feasible once sufficient computing resources 
became available. Using the approach of Ref. [40] we evaluated the transition form factors using quenched and 
dynamical Wilson quarks going to smaller quark masses than those of Ref. [40] but only at the lowest g 2 -value allowed 
on our lattices [41]. However, although we increased statistics, the quadrupolc form factors still had large errors. 
Constructing optimized sources that led to a large sample of statistically independent measurements for a given q 2 - 
value and carrying out sequential inversions through the source instead of through the current, we were finally able to 
obtain both quadrupole form factors to sufficient accuracy for a range of g 2 -values. This calculation, carried out in the 
quenched approximation, confirmed a non-zero value with the correct sign for both the quadrupole amplitudes [16]. 
Using this new methodology we extend in this work the calculation to the unquenched case. Initial unquenched results 
have been reported in Refs. [17, 18, 42]. 

Lattice calculations at the physical quark mass are currently prohibitively expensive. Recently, progress in both 
hardware performance and algorithms [43] has extended the range of accessible quark masses to lower masses, bringing 
lattice calculations into the region where they can be extrapolated using chiral perturbation theory. In this work, 
we use pions with masses as low as about 350 MeV. In a previous work [16], we found a discrepancy for the CMR 
ratio between the quenched lattice results and experiment at low values of q 2 . On the other hand, for larger q 2 - 
values, the quenched results and experiment were in quantitative agreement. Our current results for CMR show 
deviations between unquenched and quenched results at the lowest momentum transfer. In fact, the previously 
observed discrepancy is reduced by our new unquenched results. This could indicate that pion cloud effects are quite 
important at low q 2 and at light pion masses, as also discussed in the framework of chiral perturbation theory [44]. 
With pion masses in the range of 300 MeV as planned in the future one hopes to make progress in reliable chiral 
extrapolations to the physical regime. Such progress has been demonstrated recently in the extrapolation of the 
nucleon axial coupling within chiral effective theory [31, 45]. 

This paper is organized as follows: Section II gives the decomposition of the N to A matrix element on the hadronic 
level, Section III and Section IV detail the lattice analysis and outline our strategy for extracting observables. Section V 
contains our results for the transition form factors. Finally, Section VI contains our conclusions and an outlook for 
further studies we intend to perform in the field. 
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II. N TO DELTA MATRIX ELEMENT 



The electromagnetic transition matrix element can be expressed in terms of the three Sachs form factors [46], 

(A( P ',s')\ jfl \N( P ,s)) =i g) 1/2 (g^) 1/2 Mp', S ')^u( PlS ), (2) 
with the Lorentz-structure 

O a , = GMl(q 2 )K^ + g E 2(q 2 )K^ + Gc2(q 2 )K^ . (3) 
The kinematic prefactors in Euclidean space are given by 



Mi 3 m A + m N . a ,b 

^ ~ (m A + m N ) 2 + Q 2 2mN ™°v*PPP , 

E2 = _ K Mi , R Q-i (n 2, m A + m N a XjS 



j^E2 IsMl . fiO-lM2\'" AT " ,Jv o: _ a iP „ > 

K aft = (7 n +6« (Q ) 2l7 5 £ CT Aa/3P P ^ 



< 2 = -6»- V) "^™" i 75 g. ( q 2 (p + P ')» q ■ (P + P'W) , (4) 

where p(s) and p'(s') denote initial and final momenta (spins), q 2 = (p' — p) 2 , and u a (p',s') is a spin-vector in the 
Rarita-Schwinger formalism. We define £l{Q 2 ) = [{m\ + m^) 2 + Q 2 ] [{m\ — m^) 2 + Q 2 ] and Q = q, Q 4 = iq° is 
the lattice momentum transfer giving Q 2 = —q 2 . The Rarita-Schwinger spin sum for the A in Euclidean space is 
given by 



^2u a (p,s)u T (p, s) 



—ij.p + m\ 



2toa 

and the Dirac spin sum 



2p a p T .p a ^ T -p T Ja 1 
d<r r + -5—2 » 7 } o7-t7t 

3to a 3mA 3 



(5) 



^ u(p, s)u(p, s) = mN . (6) 

s 

The magnetic dipole form factor, Qmi, is the dominant form factor. The electric and Coulomb quadrupole form 
factors, Qe2, and Qc2 are subdominant. They are scalar functions depending on the momentum transfer squared. On 
the lattice we can only access space- like q 2 and therefore Q 2 > 0. The reason for using the Sachs parametrization 
in a lattice computation, as pointed out in Rcf. [40], is that the Sachs form factors do not depend strongly on the 
difference between the nucleon and the A mass. 

It is customary to quote the ratios of the electric and Coulomb quadrupole amplitudes to the magnetic dipole 
amplitude, Rem (EMR) and R,sm (CMR). In the rest frame of the A, they are given by [46, 47] 

T3 QE2{q 2 ) 

Rem 

p 151 Gc2{q 2 ) (7 s 

2m A GMi{q 2 ) 

Deformation of the nucleon- A system will produce non-zero values of Rem and Rsm- Using accurate measurements 
of the longitudinal-transverse response function at q 2 = —0.127 GcV 2 it has been shown that a spherical nucleon - A 
system is inconsistent with the experimental results [12, 48]. 

III. LATTICE MATRIX ELEMENT 

The most computationally demanding part in this evaluation is the calculation of the three-point correlation function 
given by 

(G^'" JV (t 2) ti;p , ,p;r))= cxpHp'-x 2 )cxp(+z(p'-p).x 1 )r' 3 "(^|T[ x ^(x 2 ,t 2 )^(x 1 ,t 1 )x /3 (O,0)] |fi>, 

X 2 , Xi 

(8) 
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where an initial state with the quantum numbers of the nucleon is created at time zero and the final state with the 
quantum numbers of the A is annihilated at a later time t 2 . The photon couples to one of the quarks in the nucleon 
at an intermediate time t\ producing a A. For the spin-i source, x p (x,i), and the spin-| source, x A+ (x, t), we use 
the interpolating fields 



X P{ X ) = e abc [u Ta (x) C~f 5 d b (x)] u c {x) 



(9) 
(10) 

(11) 



For large Euclidean time separations t 2 — 1\ S> 1 and t\~5>\, the time dependence and field normalization constants 
cancel in the following ratio 



X* + (x) = ^e a " C {2 [u Ta (x) C la d\x)] u c (x) + [u Ta {x) C la u\xj\ d c (x)} , 
and for the projection matrices for the Dirac indices 

r, 



1 — 1 










2 \ 


n o ) 




v0 


s 



-R CT (<2,ti,p',p,r,M) = 



G^ N (t 2 ,t uP ',p,T) 



*G NN {t 2 - tl ,P, r 4 )G^ A (t 1 , P ', r 4 )G AA (i 2 , P ', r 4 ) 



G AA (i 2 , P ',r 4 ) 



,1/2 



G AA (t 2 - ti, P ', r 4 )G JVJV (ti, P , r A )G NN (t 2 , P , r 4 ) 



t 2 -*l3>l 
*1>1 



n CT (p',p,r,^) , 

where and G AA arc the nucleon and A two-point functions given respectively by 

(G™(t,p ; r)) = 2e- i "' x r^{n|r x a (x ) i)^(o ) o)|fi) 

X 

(G AA (t, P; r)) =^e-^r^(!]|T^(x, t )xf(O 1 0)|fi) . 

X 

For the Wilson fermionic action we use the lattice conserved electromagnetic current, j^(x), given by 

j^x) = Y,Qf*f{¥{x + A)(l + ln)U^ (x)i/}f (x) - ^(z)(l - -YjUHx^ix + A)} 
/ 



(12) 



(13) 



(14) 



symmetrized on site x by taking j fl (x) — » [.^(a;) + j^(x — fi)] /2, where Q/ is the charge of a quark of flavor / and 
Kf is its hopping parameter. For domain wall fermions, we use the local four-dimensional electromagnetic current 
Qf{xp^{x)"f^ip^(x)}. This is not conserved and therefore, to relate lattice and continuum results, we need the 
renormalization constants Zy, which are known. Throughout this work we choose a frame where the A is at rest. 
For these kinematics we have q = p, p' = 0. At sufficiently large time separations, t 2 — t\ and t\, R a becomes 
independent of time yielding the desired ratio TL a . 

At the hadronic level, with the inclusion of complete sets of baryonic states and the use of Dirac and Rarita- 
Schwinger spinors the ratio of Eq. (12) leads to the relations 



n CT (0, -q ; T 4 ; fi) = iAe^pig M1 (Q 2 ) , 



(15) 



n CT (o,-q ; r fe ;j) 



M 2 (P<? S kj - PkSaj)GMl(Q 2 ) 



3 iPaPkPj 

■Z (PaSkj +PkSaj) -2 

L p 



Ge 2 (Q 2 ) 



for j — 1, 2, 3 and 



{E N -m A ) ( x 3p<rPk\ r ,,,, 



n CT (o,- q; r fc ;4) = zs( ( 5 fffe -^^)g C2 (Q 2 ) , 



(16) 



(17) 



■5 



with the kincmatical coefficients 



2 rriA + ttin 

3 4m N E N 



En + m N ' 



B = 



2mA 



-A. 



A convenient method for the evaluation of three-point functions is the sequential inversion through the sink. This 
requires fixing the hadronic state at ti to be the A with fixed vector index a. The projection matrices T at the sink 
are also fixed, but the operator inserted at any time t\ can be left arbitrary. Therefore, with one sequential inversion 
one can evaluate the three-point function for a large set of lattice momentum transfer values q, any current direction 
fi and any operator insertions at any intermediate time t\. One then looks for a plateau as a function of t\ that 
determines 11^. Eqs. (15-17) constitute a system of equations for the form factors Qmi,Qe2 and Qci at each value of 
Q 2 . It must be noted that while three independent measurements of n CT (q ; T; fi) suffice for the determination of the 
form factors, increasing the combinations of \i and photon momentum q which are measured improves the statistical 
accuracy of the form factors. With that goal in mind, we observe that Eqs. (15-17) are identically zero for several 
values of q = — p, e.g. Eq. (15) is zero when j = a or j = fi. Furthermore, for a given selection of a, not all lattice 
rotations of q are giving non-zero contributions. We therefore search for the linear combinations of Eqs. (15-17) 
which maximize the number of non-zero q contributions in a lattice rotationally invariant fashion and construct the 
following optimal combinations 



Si(q;/z) = E 11 ^ ' -q ;T 4 ;^) = iA \ (P2 -Pz)h,n + (P3 -Pi)&2,n + (Pi -P2)h tl i \Qmi{Q 2 ) 

CT=1 L > 



(18) 



3 

E 



n CT (0,-q ;T k ;fj,) 



for /j, = 1, 2, 3. For \i = 4 we have 

<S2(q;M 



-3i4|((p2 +Ps)*l,^ + (P3 +Pl)S2,n + (Pi +P2)S 3 ,„)gE2(Q 2 ) 

E N - m A 



4) 



-^\{P\P2 



P\P?, + P2P2, 



9e2{Q z ) 



2toa 



Qc2(Q z ) 



-i6B 



(PlP2 +P1P3 + P2PZ)GC2{Q 2 ) 



,(19) 



(20) 



The three-point functions involved in S\ and ^2 require one sequential inversion if one uses the appropriate linear 
combination directly in the construction of the A sink. Thus, with two inversions we obtain the maximal number of 
lattice determinations of the form factors for all the allowed lattice photon momenta q. The Si-type matrix element 
determines Qm\ while the quadrupolc form factors are extracted from the 5*2-type matrix element. The method is 
clearly superior to the method used in Refs. [40, 41] since the same CPU cost allows the evaluation of the form factors 
at all Q 2 . It should be noted that the evaluation using the SVtype sink does not determine Gc2 at the lowest allowed 
photon momentum q = (1, 0, 0) 2tt / La (and at the equivalent momentum in the other two directions). For this reason 
we use, in addition, the combination 



5 3 (q; fi) = n 3 (0, -q; r 3 ; fi) - (IT^O, -q; r l5 /1) + II ( 0, -0, -q; T 2 ; M ))/2 
to get Qc2 for the values of q for which the SVtype sink vanishes. This linear combination gives 

3 A 



<$3(q; aO 

for /z = 1, 2, 3 and 



P» 



2 

pi 



S 3 (q;/i = 4) = 



3 i B 



1 - 3 



Pi 



2mA 



Gc2(Q 2 ) ■ 



Qc2{Q 2 ) 



(21) 



(22) 



(23) 



The full set of lattice measurements for the type Si , S2 and S3 matrix elements for all contributing values of /z and 
q at a given Q 2 are analyzed simultaneously. We denote by P(q; n) the lattice measurements for the ratios R a using 
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Wilson fermions 


V 


# of confs K 


m, (GeV) 




m N (GeV) 


m A (GeV) 






yuencned, p = 


= 6.0, a 


= 2.14(6) GeV 




32 3 


x 64 


200 0.1554 


0.563(4) 


0.645(9) 


1.267(11) 


1.470(15) 


32 3 


x 64 


200 0.1558 


0.490(4) 


0.587(12) 


1.190(13) 


1.425(16) 


32 3 


x 64 


200 0.1562 


0.411(4) 


0.503(23) 


1.109(13) 


1.382(19) 






k c =0.1571 


0. 




0.938(9) 








Unquenched, (3 - 


= 5.6, a' 1 


= 2.56(10) GeV 




24 3 


x 40 


185 [51] 0.1575 


0.691(8) 


0.701(9) 


1.485(18) 


1.687(15) 


24 3 


x 40 


157 [51] 0.1580 


0.509(8) 


0.566(12) 


1.280(26) 


1.559(19) 


24 3 


x 32 


200 [52] 0.15825 


0.384(8) 


0.453(27) 


1.083(18) 


1.395(18) 






K C = 0.1585 


0. 




0.938(33) 





TABLE I: Parameters and number of gauge field configurations used in the calculation using Wilson fermions. 
the three-point function constructed using the sink-types described in Eqs. (18-21) and with Wk their statistical error. 

(Qm\ \ 
Ge2 we are led to the overcomplete set of equations 
GC2 J 

P(q;/i) =D(<i;n)-F(Q 2 ). (24) 

The kinematical prefactors, D(q, /i), are known analytically and we have inserted the continuum expressions as detailed 
in Eqs. (3) and (4). As already pointed out, we have lattice measurements for all possible vector current components 
H and photon momentum vectors q which contribute to a given value of Q 2 . For A^ such measurements, the matrix 
D(q; jj,) is an A" x 3 matrix of kinematical coefficients. The solution vector is determined from the minimization of 
the total x 2 



N /v^3 

,2 



Defining D' k - = Dkj/wk and P k = Pk/wk, (k — 1,2, ..N, j — 1,2,3) the solution is obtained using the singular value 
decomposition of the matrix D' k - (q; fi) 

D' = U-diag(X 1 ,X 2 ,X 3 )-V T , (26) 

where U is an A^ x 3 matrix, V is a 3 x 3 matrix and the Xj are the non-negative, singular values of D'. The form 
factors are therefore given by 

F = V ■ diag(l/X lt 1/A 2 , 1/A 3 ) • (U T ■ P'). (27) 

This strategy has been developed in [49] and also applied in [50] and subsequent publications. The errors Wk in the 
lattice measurements and the final error on the form factors are determined from the jackknifc procedure. Having 
decided to use sequential inversions through the sink, the sink-source separation must be kept fixed. One would like to 
use the smallest sink-source separation that guarantees that excited state contributions are negligible. In our previous 
work we found that a sink-source separation, of about 5 GeV -1 is sufficient. We check in the next Section that this 
is also sufficient for this calculation by comparing to the results obtained when we increase t 2 by about 25%. 

In the hybrid action approach we have applied hypercubic (HYP)-smearing [53] to the gauge fields. We have 
performed the computation using two different boundary conditions (b.c.) in the temporal direction. In the first case 
we impose Dirichlet b.c. on time slices and 31 using only the first half of the lattice and in the second we use the 
full lattice with antipcriodic b.c. consistent with the b.c. used in the production of the gauge fields. The lattice 
spacing a = 0.1241 fm has been determined from heavy-quark spectroscopy [54] with a statistical uncertainty of 2%. 
For Wilson fermions we use antiperiodic b.c. 

To improve convergence to the nucleon and A ground states we use Gaussian or Wuppertal smearing to create 
smeared quark fields as described in Refs. [36, 55] with smearing parameters a = 4 and n = 50. 

The parameters and number of gauge field configurations are summarized in Table I for Wilson fermions and in 
Table II for the hybrid action. 
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V 


# of confs 


Lsm-n (fm) 


(am u , d/a ) Asqtad 


(am u , d ) uwtl 


m w (GeV) 




m N (GeV) 


m A (GeV) 


20 3 x 32 


150 


2.5 


0.03/0.05 


0.0478 


0.606(2) 


0.588(7) 


1.329(9) 


1.662(21) 


20 3 x 32 


150 


2.5 


0.02/0.05 


0.0313 


0.502(4) 


0.530(11) 


1.255(19) 


1.586(36) 


28 3 x 32 


118 


3.5 


0.01/0.05 


0.0138 


0.364(1) 


0.387(7) 


1.196(25) 


1.561(41) 


20 3 x 64 


200 


2.5 


0.03/0.05 


0.0478 


0.594(1) 


0.585(7) 


1.416(20) 


1.683(22) 


20 3 x 64 


198 


2.5 


0.02/0.05 


0.0313 


0.498(3) 


0.525(8) 


1.261(17) 


1.589(35) 


20 3 x 64 


100 


2.5 


0.01/0.05 


0.0138 


0.362(5) 


0.401(13) 


1.139(25) 


1.488(71) 


28 3 x 64 


300 


3.5 


0.01/0.05 


0.0138 


0.353(2) 


0.368(8) 


1.191(19) 


1.533(27) 



TABLE II: Parameters and number of gauge field configurations used for the hybrid action. 



Domain wall fermions (DWF) [56, 57, 58] introduce an additional fifth dimension of length L 5 . They preserve the 
Ward-Takahashi identity [59] even at finite lattice spacing in the limit L5 — > 00. At sufficiently small values of the 
lattice spacing a, the effect of a finite value of L5 can be parameterized by an additional residual mass term in the 
Ward-Takahashi identity [60, 61]. This behavior describes a residual explicit breaking of chiral symmetry which can 
be minimized by choosing a sufficiently large extra dimension, L 5 . We have found that a value of L 5 — 16 in lattice 
units is sufficient to keep the residual mass, (am) roa , at most one order of magnitude smaller than the domain wall 
quark mass, (am)° WF . The height parameter of the domain wall action has been chosen to be amo = 1.7. 

The DWF quark masses displayed in Tables II have been tuned by adjusting the lightest pseudoscalar meson in the 
Asqtad calculation [27] to have the same mass as the pseudoscalar meson using domain-wall fermions. For technical 
details of this tuning procedure, see [29, 34]. 

IV. EXTRACTION OF FORM FACTORS FROM LATTICE RESULTS 

In this section we discuss several technical issues, beyond the general methodology described in the previous section, 
that must be addressed before the N to A form factors can be extracted reliably from lattice measurements. 

As mentioned already, given that sequential inversions are the most time consuming part of the calculation, the 
method of choice to calculate form factors that are functions of the momentum transfer squared is to perform sequential 
inversions through the sink. However, this approach requires that we fix the initial and final hadron states as well as 
the sink-source time separation, ti. Changing t 2 requires a new sequential inversion. 

Therefore we must first determine the optimal value of ti. The criterion is to choose ti as small as possible so 
that statistical errors due to the exponential decrease of the signal are minimized but large enough so that we ensure 
that excited states with the nucleon and A quantum numbers are suppressed. In order to decide on the value of t 2 
we compare results at two different sink-source separations in the quenched theory and in the hybrid approach. The 
quenched theory is used to provide a guideline because the statistical errors are the smallest and therefore deviations 
due to contamination from excited states can be seen more easily. We use two time separations, namely til a = 11 
and 13 and perform the calculation at the smallest pion mass. In Fig. 1 we show results for the ratio Rs x defined as 
in Eq. (12) but using the optimal source of Eq. (18) for the three-point function and normalized with the appropriate 
kimematical factors such that the plateau value yields Gmi- Note that this is not what is actually done in the 
overconstrained analysis where the plateau value for each value of the momentum vector is extracted. Nevertheless 
i?Si shown in Fig. 1 gives an idea of the quality of the plateaus that are obtained. We show results for the four 
lowest Q 2 -values when the sink-source separation is taken til a = 11. We compare the results obtained for sink-source 
separation t 2 /a = 13 at the two lowest Q 2 values and at the largest value shown in the figure. The third lowest Q 2 
value corresponding to n 2 = 3 for t 2 /a — 13 is not included since it will make the figure difficult to read. As can 
be seen the two time separations yield consistent results over a time range t\. Although in this comparison we use 
the same number of configurations the statistical errors are much smaller for the shorter time separation. Fitting to 
a constant over the plateau range but now for individual momentum vectors within the overconstrained analysis we 
obtained Gmi- The results are shown in Fig. 1, for the two different time separations using 100 configurations in each 
case. As can be seen there is very good agreement showing that a time separation of about 5 GeV -1 is sufficient. 

It is important to ensure that what we find in the quenched case carries over to the hybrid action. We expect 
pion cloud effects to be important for dynamical quarks and therefore we must check that the time evolution is large 
enough to allow the pion cloud to fully develop. Guided by our findings in the quenched theory we choose t 2 /a — 8 or 
about 5 GeV -1 and ti/a = 10 or about 6.3 GeV -1 . In Fig. 2 we show the ratio Rs 1 for the two heaviest quark masses 
whereas in Fig. 3 we show the results for Gmi at the smallest quark mass for the two different time separations. As can 
be seen, Rs 1 at the two lowest Q 2 values yields the same plateau value for both time separations. For t 2 /a = 10 the 
statistical errors are larger associated with the larger time separation. The values extracted for Gmi at the smallest 



8 



2. Or 



1.5 ■ 



m =0.41 GeV 



o 
o 



4 



$ I 5 I 



1.0 ■ 



0.5 ■ 



0.0 L 



ft =s 

o 



4 J $ 5 



2.0r 



1.5 - 



1.0 - 



0.5 ■ 



0.0 L 



0.0 



x n s =l * n 2 =3 
o n 2 =2 A n a=4 



2 3 4 5 

t, (GeV 1 ) 



* t/a = ll 

O t/a=13 



ft 



1 % ? 



0.5 



1.0 1.5 

Q 2 (GeV 2 ) 



2.0 



o 2 1 

o 

.2 
o 

4 



2 - 



m =0.60 GeV 



o 
* 



*- 



X n*=l 

o n 3 =2 



* -* 

*n a =3 

A n 2=4 



o 



* 
* 



o 

D 



m =0.50 GeV 



QfJ- -e^- -e5- -cB ^ 
*- - -*- 



* - -* 



5 
□ 



2 3 4 

t, (GeV) 



■ 



FIG. 1: Upper graph: The ratio from which Qmi is ex- 
tracted as a function of time separation for the four lowest 
Q 2 -values taking the sink-source separation t% = 11a in 
the quenched theory at pion mass m-K = 411 MeV. We 
also show results obtained with sink-source separation 
ti = 13a, at the two lowest values of Q 2 (rhombus and 
open squares respectively) and at the fourth lowest (open 
triangles), displaced to the left for clarity. The dashed 
lines are fits to the plateaus obtained with = 11a 
and span the range of fitted points. The lower graph 
shows Qmi at m w — 411 MeV, for sink-source separa- 
tions ti — 11a (asterisks) and ti = 13a (open circles), 
displaced to the right, for clarity. 



FIG. 2: The ratio from which Qmi is extracted for he 
hybrid action using a spatial volume 20 3 as function of 
time separation. We show results for the four lowest Q 2 - 
values taking the sink-source separation £2 = 8a. For the 
two lowest values of Q 2 , we show results obtained with 
sink-source separation ti — 10a(rhombus for lowest and 
open squares for second lowest momentum), displaced 
to the right, for clarity. The upper graph is for = 
0.594(1) GeV and the lower for m w = 0.498(2) GeV. The 
dashed lines are fits to the plateaus obtained with ti — 8a 
and span the range of fitted points. 



pion mass, where pion cloud effects are expected to be the largest, are also consistent for the two time separations, 
as can be seen in Fig. 3. The deviations observed at Q 2 > 1.5 are due to the large statistical noise associated with 
the larger time evolution. We therefore conclude from this analysis that a time separation of about 5 GeV -1 is also 
sufficient for our unquenched study. As a result we fix £2 ~ 5 GeV -1 or in lattice units to 11 in the quenched case, 
to 12 for Np = 2 dynamical Wilson fermions and to 8 in the hybrid approach. 

Another potential source of a systematic error is the spatial size of our lattices. Given that for the quenched case 
we use a lattice of spatial size of about 3 fm we expect finite volume effects to be negligible. Since we do not have 
dynamical Wilson configurations on a larger volume we test for finite size effects in the hybrid scheme for which, at 
the smallest quark mass, there are MILC configurations for spatial lattice size L s — 2.5 fm and L s — 3.5 fm giving 
Lgm^ = 4.6 and Lgm^ — 6.4, respectively. In Fig. 4 we show results for Qmi for these two spatial sizes. Results on the 
smaller lattice are consistent with results on the larger lattice. This indeed shows that finite volume effects are small 
for Lsm-K > 4.5. Since for all our quark masses, except the lightest mass dynamical Wilson fermions L s ra^ > 4.6, we 
expect finite volume effects to be small. 
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V. RESULTS FOR THE NUCLEON TO A TRANSITION FORM FACTORS 



The lattice results for the dominant dipole form factor are shown in Fig. 5 as we change the quark mass in the three 
types of simulations considered in this work, namely in the quenched theory denoted by Np — 0, for two degenerate 
flavors of dynamical Wilson fermions, denoted by Np = 2 and in the hybrid scheme. They are also given in Tables III, 
IV and V of the Appendix. All the results discussed in this section in the hybrid approach are obtained on the lattices 
of temporal extent 64 and using antiperiodic boundary conditions in the temporal direction. 

In all three cases the magnetic dipole decreases with the quark mass. The dashed lines are fits to an exponential 
Ansatz / exp(— Q 2 /m 2 ) whereas the dotted lines to a dipole Ansatz <7o/(l + Q 2 / m o) at the lightest quark mass. 
As can be seen, both provide a good description to the lattice results. The dipole mass that we find in the hybrid 
approach at the smallest pion mass is mo = 1.30(3) GeV. For comparison, a dipole fit to the experimental results 
yields mo — 0.78 GeV reflecting the faster fall off of the experimental results. In Fig. 6 we compare quenched and 
unqucnched results for Qmi at the three quark masses. For similar pion mass the results are in agreement even for 
pion mass as low as 350 MeV. As can be seen, lattice results fall approximately on the same curve having a weaker 
Q 2 dependence than the experimental results. In the momentum range considered here experimental results can be 
well described by a dipole form. In fact, whereas the experimental results fall off faster than the dipole form factor 
of the nucleon Gp> = 1/(1 + Q 2 /0.71) 2 with mo = 0.78 < \/0.71, the lattice results display a weaker Q 2 dependence 
yielding larger values for mo- It remains an open question whether decreasing the quark mass towards the physical 
limit will modify this Q 2 dependence. 

In Figs. 7 and 8 we show the results for the EMR and CMR ratios at the lightest quark mass, whereas we give the 
complete set of numbers in Tables III, IV and V. In the case of Wilson fermions we use both sink types S2 and S3 
to extract Ge2{Q 2 ) arid Gc'2(Q 2 )- For domain wall fermions, for which inversions are very expensive, we only use £2, 
which means that Gc'2(Q 2 ) is not determined at the lowest momentum transfer as well as at some intermediate values. 
However, we obtain results at the second lowest momentum transfer squared which for the lowest pion mass on the 
28 3 lattice is about the same as the lowest momentum on our quenched lattice, enabling a nice comparison. The main 
conclusion of this study is that these ratios are non-zero and negative. As can be seen in Fig. 7, there is agreement 
between the quenched results for EMR and those obtained within the hybrid approach at the smallest quark mass. 
For the CMR ratio there is a notable difference: At the lowest value of Q 2 the quenched results yield a value that it 
is clearly negative but smaller in magnitude than experiment. Results in the hybrid approach at the smallest quark 
mass are negative and larger in magnitude decreasing the gap between the lattice value and experiment. This confirms 
that pion cloud contributions are significant at small Q 2 modifying the quenched results. This is consistent with the 
expectation from chiral perturbation theory [44] that pion cloud contributions become significant at low Q 2 . Large 
pion contributions are also needed in phenomenological approaches [66, 67, 69]. 
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VI. CONCLUSIONS 



We have presented results for the N to A electromagnetic transition form factors calculated within full QCD, 
using two degenerate flavors of dynamical Wilson fermions, and using a hybrid action with staggered sea quarks and 
domain wall valence quarks. Comparison of these results with previous quenched calculations [16] do not reveal large 
unquenching effects in the case of the dominant dipole form factor, Qmi, for pion masses down to about 350 MeV. All 
lattice results for Qmi at these quark masses show a weaker Q 2 dependence than experiment. Comparison of results 
with the hybrid action at two different volumes indicates that volume effects are not larger than our statistical errors. 
Therefore we can not attribute this weaker Q 2 dependence to finite volume effects. Agreement between results using 
dynamical Wilson fermions and Domain Wall fermions that have different finite lattice spacing systematics suggests 
that it is unlikely to be due to discretization errors. However, we cannot presently extrapolate the lattice results to 
the continuum limit. Such an extrapolation would require three different lattice spacings at similar volumes and quark 
masses, which are currently not available. For the EMR and CMR ratios, on the other hand, unquenched lattice results 
begin to show agreement with experiment. One of the main results of this study is the quark mass dependence of the 
CMR ratio at low Q 2 . Whereas quenched lattice results underestimate the magnitude of CMR at low Q 2 , results in 
full QCD become more negative, bringing lattice results closer to experiment. This demonstrates the importance of 
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full QCD dynamics in the quadrupole form factors as one enters the chiral regime. 

The largest conceptual question as we enter the chiral regime in full QCD, is how to fully incorporate the physical 
effect of the decay of the A into a pion and nucleon on the transition form factors. Even when the pion is sufficiently 
light that the A could decay in an infinite box, the p-wave pion in a small box can still be above the decay threshold. 
Naive extrapolation to the chiral limit without including the decay threshold could produce qualitative errors in form 
factors, such as the discrepancy presently observed between the lattice and experimental magnetic form factors. In the 
event that chiral perturbation theory is sufficiently accurate at and above the A decay threshold, finite volume chiral 
perturbation theory would be an attractive framework for addressing decay channel physics. Otherwise, extension 
of finite volume techniques presently used for phase shifts will be required to address form factors. This presents an 
exciting and important research challenge. 
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VII. APPENDIX 
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-1.82(1.10) 


-3.59(93) 


0.706 


0.804 


;32) 


-0.26(99) 


-2.22(1.05) 


0.878 


0.687 


;26) 


-1.30(94) 


-5.44(1.30) 


1.042 


0.595 


;29) 


-0.95(1.24) 


-4.74(1.42) 


1.353 


0.475 


;34) 


-2.14(2.04) 


-7.01(1.22) 


1.501 


0.406 


;28) 






1.644 


0.359 


'30) 






1.783 


0.317 


; 3 i) 






1.918 


0.252 


;35) 






2.050 


0.276 


[45) 






2.178 


0.223 


'34) 






2.426 


0.172 


;s7) 







TABLE III: Quenched results for Q MU EMR and CMR 
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Wilson fermions 


Q z (GeV^) 


Qm 1 


EMR % 


CMR % 


N F = 2 Wilson, f3 = 


5.6, a' 1 = 


2.56(10) GeV 






= 691(8) 




0.447 


1.437(36) 


-0.76(43) 


-1.76(69) 


0.891 


0.989(31) 


-0.86(71) 


-3.03(1.04) 


1.304 


0.717(30) 


-1.44(1.33) 


-0.88(1.34) 


1.691 


0.509(40) 


-2.28(1.42) 


-5.75(1.95) 


2.058 


0.443(38) 


-2.08(1.48) 


-10.06(2.41) 


2.407 


0.341(37) 






3.060 


0.208(55) 










= 509(8) 




0.445 


1.210(42) 


-1.00(1.28) 


-0.93(1.42) 


0.892 


0.794(32) 


-5.26(1.87) 


-6.17(2.10) 


1.303 


0.521(32) 


-5.95(4.43) 


-3.49(3.44) 


1.685 


0.474(52) 


-5.65(3.47) 


-4.02(3.33) 


2.044 


0.296(32) 






2.384 


0.211(48) 










= 384(8) 




0.442 


1.066(43) 


-1.49(1.59) 


-1.27(1.57) 


0.893 


0.798(44) 


-1.83(2.10) 


-6.47(2.37) 


1.299 


0.589(37) 






1.671 


0.396(48) 






2.017 


0.244(47) 






2.342 


0.181(41) 







TABLE IV: Unquenched Wilson results for Q M i, EMR and CMR. 



15 





Hybrk 


1 action 




Q' 2 (GeV^) 


Qm\ 


EMR % 


CMR % 


Hyk 


irid action, 


a' 1 = 1.58 GeV 




= 


594(1) 




0.213 


1.850(44) 


-0.16(52) 




0.482 


1.434(36) 


-0.45(67) 


-3.52(1.16) 


0.738 


1.143(41) 


-0.52(1.42) 


-3.88(1.73) 


0.983 


0.954(42) 






1.218 


0.789(38) 




-6.25(2.78) 


1.445 


0.665(45) 






1.874 


0.471(68) 






2.079 


0.413(64) 






2.278 


0.363(72) 






2.472 


0.322(89) 






2.660 


0.262(145) 






2.844 


0.172(158) 








m„ = 


498(3) 




0.191 


1.557(46) 


-0.243(91) 




0.471 


1.177(38) 


-1.14(1.21) 


-1.96(1.45) 


0.735 


0.924(40) 


-0.56(2.53) 


-0.20(2.42) 


0.985 


0.783(46) 






1.224 


0.641(40) 






1.452 


0.565(47) 






1.883 


0.371(70) 






2.087 


0.381(96) 






2.284 


0.294(172) 






2.476 


0.260(240) 






2.662 


0.104(226) 






2.843 


0.146(120) 








m v = 


353(2) 




0.042(16) 


1.681(59) 


-0 79(881 




0.194(14) 


1.404(49) 


Q7( III s ! 




0.341(8) 


1.199(51) 


OQfl R9A} 


9 41 (~\ fil "l 

^.4:_Ll_L.U-L J 


0.482(9) 


1.070(48) 






0.619(8) 


0.930(47) 




-3 89(9 471 


0.751(9) 


0.813(48) 






1.005(11) 


0.723(53) 






1.127(16) 


0.660(53) 






1.246(16) 


0.623(60) 






1.362(23) 


0.581(62) 






1.475(51) 


0.518(89) 






1.586(18) 


0.525(84) 






1.695(35) 


0.518(94) 






1.906(65) 


0.392(116) 






2.009(26) 


0.328(87) 






2.111(13) 


0.298(94) 






2.209(20) 


0.275(128) 






2.306(23) 


0.230(94) 






2.402(34) 


0.222(124) 






2.497(34) 


0.073(108) 






2.682(20) 


0.087(170) 







TABLE V: Results in the hybrid approach for Qmi, EMR and CMR. For the smallest quark mass we include the errors in the 
determination of Q 2 since these are substantial for the small values of Q 2 allowed on this lattice. 



